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ABSTRACT 


This paper presents a detailed 3D model of a solid hydrogen storage tank based on metal hydride LaNi5 technology, utilizing 
COMSOL Multiphysics 6.1 software. The model takes into account the coupling of momentum, heat, mass and energy transfer 
within the LaNi5 metal hydride during hydrogen absorption. The main objective of the study is to analyze the temporal evolution 
of temperature and pressure within the tank as hydrogen is absorbed. In addition, the paper investigates the effectiveness of a 
cooling strategy involving the integration of cooling tubes into the tank configuration. This approach aims to enhance the thermal 
management of the storage system by dissipating excess heat generated during hydrogen absorption. 

Simulation results demonstrate the changes in temperature and pressure occurring within the LaNi5 metal during the process 
of hydrogen absorption. The implementation of an air-based cooling system emerges as an effective means of regulating the 
temperature of the storage tank, thus creating optimal conditions for hydrogen absorption processes. This understanding is essen- 
tial to the development of efficient thermal management solutions for solid hydrogen storage technologies. By comprehensively 
analyzing the thermal behavior of the LaNi5 metal hydride tank, this numerical study suggests that the efficient design of storage 


system is very important for rapid absorption of hydrogen. 
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1. Introduction 


Since the beginning of the 21‘ century, the entire world has shifted 
its interest towards searching for alternatives to traditional transport 
technologies due to the consequences of their impact on the climate 
[1]. All industrial sectors are aiming to overcome their dependence 
on fossil fuels, which are responsible for a substantial proportion of 
greenhouse gas emissions [2]. 

In the environmental context, hydrogen is a promising energy car- 
rier in the face of current challenges. it is a solution that can be used 
to replace fossil fuels in various fields, notably in the automotive 
sector [3], where it is used as a fuel to power a hydrogen fuel cell [4, 
5] or even internal combustion engines [6]. Thanks to its various pro- 
duction [7] and storage methods [8], it is considered a clean, renew- 
able, and highly powerful energy source. 

Currently, there are three main approaches for hydrogen storage. 
High-pressure hydrogen storage, which is realized by compressing 
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hydrogen gas at pressures ranging from 300 to 700 bars. Low-tem- 
perature liquid storage, it involves cooling hydrogen gas to as low as 
-253 °C to condensate it to a liquid state and solid-state hydrogen 
storage that requires the use of metal hydrides to absorb hydrogen 
gas [8]. 

The use of metal hydride-based hydrogen storage offers a signifi- 
cant volume density (of the order of 120 kg/m?) compared with the 
other two storage methods, liquid (71 kg/m’), and gaseous (42 kg/m’) 
[9]. However, the enormous amount of heat emitted during the hy- 
driding process is a major drawback [10]. This poses a challenge to 
the thermal management of the storage tank [11]. 

The storage of hydrogen in solid form utilizing metal hydrides ne- 
cessitates thorough consideration of the thermal challenges. In order 
to be enable to absorb the hydrogen the tank must be maintained at 
an appropriate temperature [12]. 

The main objective of this work is to study the absorption phe- 


Bousseksou and Lin 


Engineering Perspective 4 (1): 32-39, 2024 


nomenon of a metal hydride-based solid-state hydrogen storage sys- 
tem for powering a PEMFC (proton exchange membrane fuel cell) 
fuel cell used in vehicles. 

To optimize the performance of metal hydride tanks, research into 
tank design and configuration is essential for efficient operation of 
stored hydrogen distribution. The present work implies a numerical 
study on the comparison of different techniques proposed to improve 
the thermal management of a solid-state hydrogen storage tank 
based on a LaNi5 metal hydride. In particular, the evolution of the 
temperature within hydride bed upon absorption. 


2. Numerical model 


This study, a numerical model was used to simulate the thermal 
behavior of a solid-state hydrogen storage tank. The physical model 
of the metal hydride tank studied is composed of a solide phase 
(metal: LaNi5) and a pure gaseous phase (hydrogen: Hz), thus form- 
ing a porous medium as showed in Figure 1. The choice of tank ge- 
ometry is based on a 3D type is seen on Figure 2. The tank consists 
of a cylindrical of radius R = 0.25 m and height H = 1 m. The hydro- 
gen is injected into the tank from the top at Py2 in= 8 bars during the 
absorption of hydrogen. The lateral surface experiences cooling 
through natural convection, while the remaining orthogonal surfaces 
are presumed to exhibit adiabatic behaviour. The temperature is 
maintained constant at the extremities of the tank and set to the value 
of To = 293 K. 


Figure 1. Schematic section view of the hydrogen absorption in a porous 
media 
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Figure 2. The geometry of the model 


To simplify the numerical model, the following assumptions are 
considered [13, 14, 15, 16, 17]: 


The gas phase (hydrogen) is pure. 

The tank is considered in three dimensions (3D). 

The media are in local thermal equilibrium between gas and 
solid. 


iv. The validity of Darcy’s law. ?z 
v. The radiative transfer in the porous medium is neglected. 
Vi. The hydrogen inlet temperature and pressure are maintained 


constant (To and Po). 


The equations governing the hydriding process are discussed in 
detail below. 


2.1 Energy conservation 


According to the [14, 18] heat transfer modeling can be simplified 
by setting a temperature variable for the system (Eq. (1)). 


OT 
(PCy ere 5, + PgCpg¥g VT = V(Aere VT) + sm (1) 


Where the effective heat capacity, (PC) ) or, is given by Eq. (2): 


(PCy ert = €Pglyg + C= €)psCys (2) 
The effective thermal conductivity can be written as: 

Aepp = €Ag + (1 —€)A, (3) 
The heat source term S,, is equal to: 

Sm = (1 — €). 1 AH 11.222 (4) 


: F é ; . Ops . 
With, € is the porosity of the medium studied and es is the amount 


of hydrogen absorbed by the metal as a function of time. 
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In the case of absorption [19]: 


=Ea P 
ee Ca. err .In (=) (Dss > Ps) (5) 
Peq 
With Pq is the equilibrium pressure. 
2.2, Mass conservation 


For gaseous media (hydrogen): 


(6) 


a : 
€ m + div(py2Uq2) = —m 


Where, Py is the density of gas (hydrogen) and v4, is the veloc- 


ity field of the gas (hydrogen). 


For solid media (metal hydride): 


(1- €) = =r 


(7) 
Where, ps the density of solide (bed) and m the mass transferred 


from the gas phase to the solid in function of time. 


The gas density py, is described by the equation: 


(1— ) P= 1A 


(8) 


2.3 Momentum conservation 


According to Darcy’s law the hydrogen gas velocity within the 
thank is: 


(9) 


Where, k is the intrinsic permeability of the porous medium (m’), 
and Ln? is the dynamic viscosity of the fluid (Pa.s). 


2.4 Kinetics reaction 


The absorption kinetic of a hydride metal is related to the 
amount of hydrogen absorbed by the metal over time can be ex- 
pressed as: 


Ops 


ot 


= Cx exp (—F2) n(#) (sae ~ 6s) (10) 


The equilibrium pressure P., is calculated using the Van't Hof 
relation via the following term [19]: 


Lf 


In (2) =A- 


Po 


(1) 


With A and B are Van't Hoff constants having the values of 
12.95 and 3731.42 respectively [19]. 
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2.5 Initial and boundary conditions 


Initially "t = to", the temperature, pressure and hydride density 
are assumed to be constant (Figure 2). 


T(to, 7,2) =T) (12) 

P(to,7,Z) = Py (13) 

p(to, r, Zz) = Po (14) 
At the hydrogen inlet (Figure 2): 

T(x, y,t) =Tp (15) 

P(x,y,t) = Po (16) 
The lateral cooling wall (Figure 2): 

Reonv (T a Ty) (17) 


2.6 Modeling parameters 


The thermophysical proprieties of the LaNi5 and Hz used in this 
study are succinctly presented in Table 1. 


Table 1. Thermophysical proprieties of the metal hydride (LaNi5), hy- 
drogen (Hz) and other parameters used in the simulations [13, 14, 16, 17, 
20, 21). 


Pin 8 [bar] 

Tin 293 [K] 
M2 2.01588 [g/mol] 
An2 0.24 [W/m. K] 
pH2 0.0838 [kg/m?>] 
Cpr2 14.890 [J/mol. K] 

Po 1 [bar] 

To 293 [K] 
Heonv 1650 [W/ (m?.K)] 
Knya 2.4 [W/ (m.K]] 

€ 0.5 
Cps 419 [J/kg. K] 

k 3x10"? [m?] 
Psat 8520 [kg/m?] 

ps 8400 [kg/m?] 
CPair 1005 [J/kg. K] 
kair 0.025 [W/(m?.K)] 

Cpr0 4180 [J/kg. K] 
kr0 0.598 [W/(m’. K)] 

Tc 293 [K] 

R 0.02 [m] 

Ea 21170 [J/mol] 

Ca 59.187 [1/s] 
AH 30 [kJ/kg] 

Rg 8.314 [J/(mol. K)] 
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2.7 Cooling strategies 


In order to enhance the cooling efficiency of the storage tank, a 
novel approach was contemplated. This entailed the incorporation of 
internal cooling tubes withing the metal hydride, as illustrated in Fig- 
ure 3 (b and c). While maintaining the original geometry of the tank, 
two cooling tubes with a radius of 0.02 m were integrated. 


For this process, two heat transfer fluids were used: air and water. 
These fluids were chosen because of their specific heat transfer prop- 
erties, abundance, and their common use in cooling applications. 
The aim of this study is to evaluate the effectiveness of the afore- 
mentioned approach using two distinct fluids, compared with the 
case of a basic tank which cools solely by natural convection. A 
schematic of the different configurations is presented in Figure 3. 
The cooling tubes, which are supposed isotherm (Figure 3): 


T(t,7,Z) =Tp (18) 


LaNiS + H, 


Figure 3. Schematic section view of the tank equipped with two cooling 
tubes. (a) basic tank, (b) tank with two air-filled cooling tubes and (c) 
tank with two water-filled cooling tubes 
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3. Results and discussion 


In this section, we determine the variation of temperature and 
pressure as a function of time during the absorption process in hy- 
drogen storage tanks. Then, we will provide a comparative assess- 
ment of the temperature variation in the tank for the various pro- 
posed cooling techniques. 


3.1 Temperature evolution 


Figure 4 shows the temperature variation as a function of time at 
a specific measuring point within the tank. A rapid increase in tem- 
perature can be observed, reaching the maximum value (345 K), at- 
tributed to the exothermic nature of as the hydriding reaction of the 
intermetallic. Then, after this rapid increase, the temperature gradu- 
ally decreases towards ambient temperature (to 298 K). Our results 
align with those reported in the literature [12]. 
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Figure 4. Temperature variation as a function of time during absorption 


The spatio-temporal evolution of the temperature inside the thank 
from 0 to 4000 is shown in the Figure 5. The temperature rises rap- 
idly, reaching a maximum value of 345 K, followed by a gradual 
descent to reach a peak temperature of 330 K after 4000 s. 
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Figure 5. Spatio-temporal temperature variation 


3.2 Pressure evolution 


Figure 6 illustrates the temporal variation of hydrogen pressure 
within the hybrid bed. It is observed that the pressure increases 
rapidly until it reaches a steady state corresponding to the set pres- 
sure. This rapid and uniform stabilization of pressure is observed 
throughout the entire hybrid zone, leading to the conclusion that 
H2 percolation times are extremely short [22]. 


Measurement} 
point 


Pressure (bars) 


8000 


1000 2000 3000 4000 5000 6000 7000 


Time (s) 
Figure 6. Temperature Pressure variation as a function of time 


3.3 Effect of the tubes cooling 


Figure 7 illustrated the temperature field reigning inside the tank 
from 0 s to 4000 s in three different configurations ((a) basic tank, 
(b) tank with two air-filled cooling tubes and (c) tank with two water- 
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filled cooling tubes). We can observe the rapid growth of the tem- 
perature in the hydride bed at the start for all three tank configura- 
tions. Then, they decrease progressively as the kinetics of the hydrid- 
ing reaction slow. However, the rate in which the temperature goes 
down shows a marked disparity between configurations with cooling 
tubes and those without, where cooling is based solely on the periph- 
eral wall of the tank. We observe a rapid increase, reaching a maxi- 
mum value of 345 K for t= 500 s, both for the configuration without 
cooling tubes and for that with air-filled tubes. However, in the case 
of the configuration with water-filled tubes, a slightly lower maxi- 
mum temperature of 344 K is observed. Thereafter, it gradually de- 
creases to reach a temperature of 330 K after 4000 s in the case of 
the tank without cooling tubes, while for the configurations with 
cooling tubes filled with air and water, the corresponding tempera- 
tures after 4000 s are 319 K and 318 K, respectively. 

The presence of cooling tubes in the storage tank reduces the time 
required to cool the hydride bed, which in turn reduces the hydriding 
time. This approach considerably increases the heat exchange sur- 
face area within the tank. The cooling tubes also enable the heat to 
be dispersed more efficiently. 

The use of water as a coolant is highly effective due to its excel- 
lent heat transfer properties, making it an efficient cooling medium. 
Nonetheless meticulous selection of materials for the design of cool- 
ing tubes is imperative to mitigate potential issues such as corrosion 
and leaks, thereby ensuring the enduring robustness of the envisaged 
cooling system. However, it's worth noting that in some cases, air 
may be preferable to water as a cooling medium. While water is 
highly efficient for heat transfer, it can also present additional chal- 
lenges in terms of corrosion and maintenance, especially if the con- 
struction materials are not compatible with water. Conversely, uti- 
lizing air as a coolant offers a more straightforward and potentially 
more reliable alternative, exhibiting reduced susceptibility to opera- 
tional issues. This can be particularly advantageous in applications 
where water-related risks are high or where simplicity and ease of 
maintenance are priorities. 

Therefore, considering that there is only a 1 K temperature dif- 
ference between air cooling and liquid cooling, opting for air cooling 
may be more effective compared to the potential issues associated 
with using water. 
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Figure 7. Variation of the temperature in function of time during absorption for different configurations: (a) basic tank, (b) tank with two air-filled cool- 
ing tubes and (c) tank with two water-filled cooling tubes 
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4. Conclusions 


This study contributes to a better understanding of hydrogen 
storage in tanks, in particular the hydrogen absorption process 
with different designs of the cooling system. A three-dimensional 
mathematical model describes how heat and mass transfer in 
metal hydride are carried out. The results showed a substantial and 
swift generation of heat during the hydrogen absorption. Further- 
more, this study compares cooling solutions for solid hydrogen 
storage tanks. It offers different approaches to controlling the tem- 
perature and thus ensuring safe hydrogen storage. Incorporating 
cooling tubes filled with air significantly enhances heat transfer 
within the tank. 


There are many perspectives and avenues for achieving our 
study: 


i. Desorption phenomenon: In the current work we have 
studied the phenomenon of hydrogen adsorption. Subse- 
quently, a systematic examination of hydrogen desorp- 
tion kinetics is envisaged for future research. This will 
involve examining the influence of key parameters in- 
cluding temperature, pressure, material characteristics 
and surface morphology on the desorption process (a 
Model utilizing COMSOL Multiphysics is currently un- 
der development). 


Incorporating phase-change materials (PCMs) within hydrogen 
tank facilitates the storage of hydrogen at elevated densities and 
reduced temperatures when compared to high-pressure and cryo- 
genic tanks. The integration of a PCM is an important part of the 
thermal management of the tanks and offers several important 
advantages, such as the recovery of excess heat and the mainte- 
nance of temperature at optimum levels. There are many per- 
spectives and avenues for achieving our study. 
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Nomenclature 
Pin Inlet pressure of hydrogen 
Tin Inlet temperature of hydrogen 
M2 Molecular mass of hydrogen 
An2 Heat conductivity of the Hz gas 
PH2 Density of the H2 gas 
Cpre2 Specific heat of H2 gas 
Po Initial pressure 
To Initial temperature 
heonv Heat convection coefficient 
Knya Specific heat of the solid 
€ Permeability of the metal 
Cps Specific heat of the solide 
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k Permeability of the metal 

Psat Saturated metal hydride density 
Ds Ho-free metal hydride density 
CPair Specific heat of air 

kair Heat conductivity of air 

Cpi20 Specific heat of H20 

ki20 Heat conductivity of air 

Te Cooling tubes temperature 

R Annular disc unit radius 

Ea Activation energy for absorption 
Cs Absorption rate constant 

AH Reaction heat of formation 

Re Universal gas constant 
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